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Abstract. We consider a smooth counting function of the scaled zeros of 
the Riemann zeta function, around height T. We show that the first few 
moments tend to the Gaussian moments, with the exact number depending on 
the statistic considered. 



1. Introduction 

In this paper we will examine linear statistics of zeros of the Riemann zeta 
function. Denote its nontrivial zeros by 1/2 + rjj, j = ±1, ±2, . . . with 7_j = — 7j 
and !SHe(7i) < 0^(72) < .... Let N(T) denote the number of zeros in the strip 
< $He( 7 ) < T, then N(T) = N{T) + S{T) where 



N(T) = 1 + -3m log U- iT / 2 Y(\ + \iT) 



7T 

T T 7 
= — log— + -+0(1/T) 
2n zire 8 

Selberg fj has studied the distribution of the remainder term in the counting 
function, S(t) = N(t) ~ N(t), as t varies between T and T + H, where H = T a 
with 1/2 < a < 1. He showed that S(t) has Gaussian moments in the sense that 
when T — ► 00, 

2k 



1 

H 



T+H 



S(t) 



„ (2*)! 



v/OogiogW^ 

Fujii 0], among others, has studied the distribution of N(t + h) — N(t) around t 
near T. Since y L S(t) dt — > 0, the mean of this is asymptotic to N(T+h) — N(T), 
and the error term, S(t + h) — S(t), (which is thus asymptotically centered) has 
Gaussian moments, so long as h is larger than the mean spacing of zeros at that 
height. That is, if h log T — > 00 but h -C 1 then 

where 

/ll °s T 1-cost , 

dt 

t 

He has similar results for when h — > 00 subject to /i -C T, when cr 2 is replaced by 

cr 2 = 1 (io g iogr-iog|c(i + ift)|). 
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Note that if h is of the order of the mean spacing, that is if h = 0(1/ log T), 
then the main term is the same size as the error term (that is, both are 0(1)), and 
we may no longer conclude the distribution is Gaussian. This is not surprising, 
since for h = 0(1/ log T) the distribution of N(t + h)- N(t) in the large T limit is 
discrete. 

In this paper we will study the counting function in that critical scaling. Rather 
than study N(t) itself, instead we will investigate the distribution of a smooth 
version of the counting function in intervals of size comparable to the mean spacing, 
2-71"/ logT. In particular, for a real- valued even function /, and real numbers r and 
T > 1, set 

Nf(r):= /(^(7i-T)). 

j=±l,±,2,... 

If / is the characteristic function of an interval [—1, 1] and if all the jj are real, 
then Nf (r) counts the number of zeros in the interval [r — 2ir/ log T,t + 2tt/ log T] . 
However, we will take / so that its Fourier transform, f(u) := /_ f(x)e~ 2mxu dx, 
is smooth and of compact support, and will not assume the Riemann Hypothesis. 

As T — > oo, we consider the fluctuations of Nf(r) as r varies near T in an interval 
of size about H = T a , where < a < 1. More precisely, given a weight function 
w > 0, with J_ w(x)dx — 1, and w compactly supported, we define an averaging 
operator 

We will show that the expected value of Nf is 

J —OO iU & A 

We will also show that for / e C C °°(M) the first few moments of ((N f ) m ) TH of 
Nf are Gaussian: 

Theorem 1.1. Let H = T a with < a < 1, and let f S C^°(E) be such that 
supp / C (— 2a/m, 2a/m). Then the first m moments of Nf converge as T — * oo to 
those of a Gaussian random variable with expectation f(x) Ax and variance 

2 



(1) aj= / min(|tt|,l)/(u) J du. 

J — oo 

The local statistics of the critically scaled zeros of the Riemann zeta function 
around height T (that is, zeros scaled by the mean density, ^r-) are believed [Q, || 
to behave like eigenangles of a random unitary matrix, when scaled by N/2ir, which 
is their mean density. Indeed, a similar result to the theorem above holds in random 
matrix theory Since the 6 n are angles, we consider the 27r-periodic function 



oo 



F N (9):= J2 f(U° + ^)) 

3- 



and model Nf by 



JV 

Z f (U) :=^Fjv( 



where U is an N x N unitary matrix with eigenangles 



in- 
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Writing E to denote the average over the unitary group with Haar measure, 
then without any restrictions on the support of the function /, we found in 0] that 
E{Zf(U)} = J_ f(x) dx, and that the variance is 

/>oo 

(2) a{ff= / min(l, |u|)/(u) 2 du 



Observe that this is in complete agreement with the mean and variance of Nf if 
/ has the same support restrictions. Furthermore, we showed that for any integer 
m > 2, if supp/ C [— 2/m, 2/m], then 

x m I if m odd 

lim E{ (Z f -E{Z f } m }= < m , ., 
N ^°° \ 2^(1/2)^ ) if m even 

where <r(/) 2 , the variance, is given in (||). 

These are the moments of a normal random variable. However, the higher mo- 
ments are not Gaussian, and we called this "mock-Gaussian behaviour" 



The random matrix results suggests that theorem 1.1 is not the complete truth. 
We expect the variance of Nf , (|l]) , to hold without any restriction on the support of 
/, and the m-th moment of Nf to be Gaussian so long as supp / C [—2/m, 2/m]. 

Random unitary matrices can also be used to model the low-lying zeros of Dirich- 
let L-functions, and mock-Gaussian behaviour was found there too 0. Other clas- 
sical groups (symplectic, special orthogonal) are believed to model different classes 
of L-functions (like the quadratic L-functions), and they also show mock-Gaussian 
behaviour [Bfl. 



2. Proofs 

2.1. The explicit formula. Set Q(r) = |*(± + ±ir) + \^{\ - |ir) - log tt, where 
i]/(s) = 5r(s) is the polygamma function. We need the following version of Rie- 
mann's explicit formula: 

Lemma 2.1. Let g(u) E C c °°(i?) and let h(r) = g(u)e iru du. Then 
(3) 

• * -i /* OO OO A / \ 

E ^i) = + K\) + ^ j ^ HrMr) dr-E^ i^gn) + g(- logn)) 

where A(n) is the von Mangoldt junction. 

For / E C C °°(IR), setting 

loeT e _i ™ - u 

Kr)=f{-l-{r-T)), g( u ) = ——/(- > 



2tt v av ' logT J UogT 7 

we have Nf(r) — Nf(r) + S/(t) where 

+ /r^-,)V/f^(4-) 



2vr 2 '/ V 2tt 2 
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and 

(5) Sf{T) = -^g^/^) ( e — + e — ) 

Remark. The conditions on / (which are determined by the explicit formula), that 
its Fourier transform has compact support and is infinitely differentiable, can be 
considerably weakened to requiring that f(r) is analytic in the strip — c < 3m(r) < 
1 + c (c> 0) such that f(r) < (1 + \r\)~( 1+ V (S > 0, r e M). 

2.2. The mean. 

Lemma 2.2. For aZZ A > 1, 



[°° f( X )dx+o( 1 W ) +0 

'-oo 1 + |r| (logT) 2 



logTy_ 00 ' w v l + |r| (logT) 2 ' V (logT)^ 4 

Proof. To evaluate the integral, we change variables and split the domain of inte- 
gration into two parts: 

, 2tt V " 1 ; 2tt logT J ; V logT^ 



^(/ +/ + -^)dx 

g^ '|x|<Y '|x|>Y lo g T 



logT J\x\<Y J\x\>Y 

where Y -> oo but Y = o(logT), say Y = VTogT. 

For the bulk of the integral xj logT is small, and so we expand 

to find 



logT J\ X \<Y 

= logT J_ 00 Jy '"^ ■ ~ y l + \r\ (logT)S 

For the tail of the integral we use /(x) <C |x|~ w for any N ^> 1 (which follows 
from / e C C °°(R)) and Stirling's formula (which yields fi(r) = log(l + |r|) + 0(1) 
for all r £ E) to find that it is dominated by 

2JV , ,„ Y 2./V /" 1 1 



FJV-ilogT fav 11 logT' \ogTJ x>Y x"- 1 (x+(l + |r|)logT) 

„ log(l + |r|) 

yw-i logT 
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Thus we have 



f°° logT dr n(r) f°° 

/ f{~^ — (r-T )mr)— = - — -/ /(x)da 



4+|r| (logT) 2j V^-ilogT^ 
Taking F = ^/logT gives (|). □ 

Lemma 2.3. For all f £ C C °°(M), 

(AT/) = / /(x)d.x + T - oo 

Proof. Since fi(r) = log(l + |r|) + 0(1) for all r, we have 

fifr) Z" 00 \ f°° /(x)da; / f°° t-T dr 



logT7_oo / T)jr logT Vi-oo * H 'H 

logT 



and since the average of the error term in lemma 2.2 is similarly 0(1/ logT), we 
have that 

s £ ' - ->) °m*)„ - £ + °<i^' 

The averages of the polar terms, /(^(^ - r)) + /(^(-j - t)), in JV/ are 
bounded by 0( gl ^ gT ) since by Parseval 

, M 2tt v 2 " y H 'H J_ OQ \ogT JV logT' v w; y 

and since u; has compact support, the integral is over |y| <C 1/H and is bounded 
by 0{l/H\ogT). ' ' n 

Proposition 2.4. For / with f £ (7^°(M), if H oo then the mean value of Nf 
is given by 

(7) {N f ) TjH = f°° f(x)dx + 0(-±-), T^oo. 

J — oo ^og 1 



Proof. In view of Lemma 2.3, it suffices to show that the mean value of Sf is zero 



as H — > oo. Indeed, we have 

W« = i^i:^0 ( S (|.o g n )e — ■ + S (-| 1 o g , l)e -. 

Since w has compact support, and the prime powers n are at least 2, the summands 
vanish once H is larger than a certain constant which depends upon the support of 
w. □ 
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2.3. The centered moments. 



Proof of Theorem 1.1. From lemma 2.2 and proposition 2.4 we have 

m / \ 

((*/ - ^ H ) m ) TH - e (:) (^-(v>^ 

n=0 v 7 



and so it is sufficient to show that the mth moment of Sf is the same as that of 
a centered normal random variable with variance given by (Q). This is achieved in 
the theorem □ 

Before we calculate the m-th moment of Sf, as a warm-up we will consider the 
variance. 

Proposition 2.5. If H = T a with < a < 1 and supp/ C (—a, a) then we have 

t°° * 1 

(( s f) 2 ) T . H = / min (M' i)/(«) 2 d« + o{— - ) 

J — CO ^& ^ 

Proof. Using the expression (||), multiplying out (Sf) 2 and integrating we find 

° £l,£2=±lll,l2 V V <= o 

x (eilogm + e 2 logn 2 ))e~ iT(eilog ™ 1+£2los " 2) 

27T 

In order to get a nonzero contribution we need e\ = — e 2 (since once H is larger 
than a certain constant which depends upon the support of w all the ei = e 2 terms 
vanish). Furthermore, since supp/ C (—a, a) we have n\ < T a ~ e for some e > 0, 
and therefore 

w(— log— ) = 
27r n 2 

(for large enough T) unless ni = n 2 . 

Therefore, taking into account that w(0) = J"_ w(a;)dcc = 1 we find as soon as 
T is sufficiently large, 



A(n) 2 j- logn . 



2 



log 2 T^ n TogT" 
We note that by the Prime Number Theorem, as T — * oo 



(logT) 



1 

-ioir ; 



where a 2 is given by (0), this being true since supp / C (—1, 1) by assumption of 
the theorem. □ 
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Theorem 2.6. If H = T a with < a < 1 and supp / C (—a, a) then for 

2 < m < 2a/a 

we have 

( (2fe)! 2fe _ 2j L i 

;;;:^ +1 +^> 

where the variance a'j is given by ([!]). 

Proof. Using the expression (|5|), multiplying out (Sf) m and integrating we find 

<w>,, -c^r E E fi^O 

6 ei,...,e m =±lni,...,n m j=l V J 6 



Since u> has compact support, in order to get a nonzero contribution we need 



| ^ e 3 - log njl < — 

Set M = lie =+i n j an< ^ ^ = Lie =-i n j- ^ 7^ then assume w.l.o.g. that 
M > N, say M = N + u with u > 1. Thus for a non-zero contribution we need 

1 , M , , u , 1 
— » log — = log(l H ) » — 



and hence T a = H ^ N < y/MN < T m ( Q ~ £ )/ 2 since ri,- < T Q ~ e by assumption on 
the support of /. But a < 2a/m, so this is a contradiction. Therefore M = N, and 
E e i " ; = 0. 

Thus for T ^ 1, we find (taking into account that w(0) = J_ w(x)dx = 1) 



to ei,...,e m =±l n 1 ,...,n m >2 j=l v to 

ETLi e j logrij=0 



(8) = E J (^) 

£C{l,..,m} 

where 

w w-(^r e n^O 

b m,...,n m >2 i=l V 3 b 



B J 



and the subset of indices E denotes the ej which are positive. That is j S E iff 

Fix a subset £c {1, . . . , m}. The sum in J(E) is over tuples (m, . . . , n rn ) which 
satisfy Jlje-E n J = l~L^£: n 2- We sa y that there is a perfect matching of terms if 
there is a bijection a of E onto its complement E c in {1, . . . , m} so that = n a fj\, 
for all j 6 J5. This can happen only if m — 2k is even and #E = #£' c = k. 

Decompose 

(10) J{E) = J diag (E) + J non (E) 
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where Jdiag (E) is the sum of matching terms - the diagonal part of the sum (nonex- 
istent for most E), and J n on(E) is the sum over the remaining, nonmatching, terms. 

2.3.1. Diagonal terms. Assume that m = 2k is even. There are ( 2 fe fc ) subsets E C 
{1, . . . 2k} of cardinality k = m/2, and for each such subset E, Jdiag(E) is the sum 
over all fc! bijections a : E — > E c of E onto its complement, of terms 

^(logT) 2 ^ n J[ \ogT ) J 
We evaluate each factor by using the Prime Number Theorem: 

1 v A(n) 2 ^ logn 2 1 f°° \ogt ? \ogt 2 

(logTpi- n n \ogT> ~ (logT) 2 7 2 t ^logT j 

uf(u) 2 du 

Since our function is even and supported inside (—a, a) and a < 2a /m < 1, we can 
rewrite this as 

i r 00 

- min(l,M)/( U ) 2 d U - :( 7(/) 2 /2 

This shows that for m — 2k even we have as T — > oo that 

)! 

2 fe fc! 



BC{l,...,m} 

Below we will show that the nondiagonal terms J n on(E) are negligible, and hence 
by (|8|) and ( |l0| ) we will have thus proved Theorem 2.G. □ 

2.3.2. Bounding the off-diagonal terms J non (E). We will show that 
Lemma 2.7. 

Jnon(E) <C 



logp 1 



logT 

Proof. Since 

logT 2^2^ p k/2 ^ \ ogT 2^ p 3/2 ^ logT 

to p fe>3 ^ P 

the contribution of cubes and higher prime powers to (|^) is negligible, and we 
may assume in J non (E) that the rij are either prime or squares of primes (upto a 
remainder of 0(1/ logT)). By the Fundamental Theorem of Arithmetic, an equality 
Llje-E n i = IIie.Ec n i forces some of the terms to match, and unless there is a 
perfect matching of all terms, the remaining integers satisfy equalities of the form 
n\U2 — TI3 with ni = n2 = p prime and 713 = p 2 a square of that prime. Thus upto 
a remainder of 0(1/ logT), J non (T) is a sum of terms of the form 



1 \ - (logp) 2 f Aogp k ]2 I [ 1 \p (logp) 3 f Aogp ^fAogP 2 - 



(l 0gT )2 Z_. p k JV logT' | \^(logT)3^ p 2 J UogT y J MogT 



\ fe=l,2 

with 2u + 3w = rn, and u > 1 
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We showed ( |11| ) that the matching terms have an asymptotic value, hence are 
bounded. We bound the second type of term by 

1 ^ (logp) 3 - logp 2 c- logp 2 1 ^ (logp) 3 1 



(logT)3 ^ v 2 log T log T (logT)3 ^ P 2 - (logT)3 

Thus as long as v > 1 (that is if there is no perfect matching of all terms), we get 
that the contribution is O ( 1 / log T) . □ 
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